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ABSTRACT 



We provide an in-depth analysis of the He II and H I absorption in the inter- 
galactic medium (IGM) at redshifts z = 2.3-2.9 toward HE 2347-4342, using spec- 
tra from the Far Ultraviolet Spectroscopic Explorer (FUSE) and the Ultraviolet- 
Visual Echelle Spectrograph (UVES) on the VLT telescope. Following up on our 
earlier study (Kriss et al. 2001, Science, 293, 1112), we focus here on two major 
topics: (1) small-scale variability (Az 10~^) in the ratio rj = N(He II)/N(H I); 
and (2) an observed correlation of high-rj absorbers (soft radiation fields) with 
voids in the (H I) Lya distribution. These effects may reflect fluctuations in the 
ionizing sources on scales of 1 Mpc, together with radiative transfer through a 
filamentary IGM whose opacity variations control the penetration of 1-5 ryd radi- 
ation over 30-40 Mpc distances. Owing to photon statistics and backgrounds, we 
can measure optical depths over the ranges 0.1 < rneii < 2.3 and 0.02 < thi < 3.9, 
and reliably determine values oi rj 4rHeii/THi over the range 0.1 to 460. Val- 
ues f] = 20 — 200 are consistent with models of photoionization by quasars with 
observed spectral indices as = — 3. Values rj > 200 may require additional con- 
tributions from starburst galaxies, heavily filtered quasar radiation, or density 
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variations. Regions with rj < 30 may indicate the presence of local hard sources. 
We find that r) is higher in "void" regions, where H I is weak or undetected and 
~80% of the path length has 7] > 100. These voids may be ionized by local soft 
sources (dwarf starbursts) or by QSO radiation softened by escape from the AGN 
cores or transfer through the "cosmic web" . The apparent differences in ionizing 
spectra may help to explain the 1.45 Gyr lag between the reionization epochs of 
H I (^Hi ~ 6.2 ± 0.2) and He II (^Heii ~ 2.8 ± 0.2). 

Subject headings: intergalactic medium — quasars: absorption lines — ultravio- 
let: galaxies 

1. INTRODUCTION 

In the last decade, studies of the intergalactic medium (IGM) have been enhanced by 
high-resolution spectra of the Lya absorption lines of H I and He II. This pioneering work 
includes spectrographic studies of H I with the Keck telescope (Hu et al. 1995) and VLT (Kim, 
Cristiani, & D'Odorico 2001) and of He II by the Hubble Space Telescope (HST) (Jakobsen 
et al. 1994; Hogan, Anderson, & Rugers 1997; Heap et al. 2000; Smette et al. 2002) and 
the Hopkins Ultraviolet Telescope (HUT) (Davidsen, Kriss, & Zheng 1996). With FUSE^ we 
are now able to obtain moderate-resolution spectra of He II absorption (Kriss et al. 2001, 
hereafter denoted KSOl). Taken as a whole, these studies confirm the general theoretical 
picture (Gen et al. 1994; Zhang et al. 1995; Miralda-Escude et al. 1996; Dave et al. 2001) 
of the IGM as a fluctuating distribution of baryons, organized by the gravitational forces 
of dark matter and photoionized by high-redshift quasars and starburst galaxies (Haardt & 
Madau 1996; Fardal, Giroux, & ShuU 1998, hereafter denoted FGS98; Madau, Haardt, & 
Rees 1999). 

Studies of absorption in the H I and He II Lya fines at 1215.670 A and 303.782 A 
respectively, are particularly important for understanding the transition of the high-redshift 
IGM from a neutral to ionized medium. In practice, the He II absorption becomes detectable 
in the far ultraviolet at redshifts ^ > 2.1 with FUSE, and at z > 2.8 with HST. The 
"reionization" epochs appear to occur at 2; ~ 6.2 ±0.2 for H I (Becker et al. 2001; Djorgovski 
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et al. 2001; Fan et al. 2002) and at z fa 2.8 ±0.2 in He II (Reimers et al. 1997; KSOl; Smette 
et al. 2002). 

The He II Lya absorption is generally much stronger than H I Lya, by a factor denoted 
by the parameter rj =N(He II)/N(H I), predicted theoretically (FGS98) to be 50-100 for 
photoionization by quasar backgrounds. For optically thin lines, we can write r) 4THeii/THi. 
The greater strength of He II arises because it is harder to photoionizc than H I, owing to 
lower fluxes and cross sections at its ionizing threshold {huT = 54.4 cV), and because He III 
recombines 5-6 times faster than H II. Variations in rj can therefore be used to diagnose the 
sources and fluctuations of the metagalactic ionizing spectrum in both the continua of H I 
(1 ryd) and He II (4 ryd). The He II absorption is also a good diagnostic of absorption from 
low-density regions ("voids") in the IGM. Because rj ^ 1 and rneii = (?7/4)rHi, He II can be 
tracked into much lower density regions than H I. In principle, He II absorption can be used 
to probe hotter, coUisionally ionized gas at T 10^ K, produced by shock-heating (Cen & 
Ostriker 1999; Dave et al. 2001). 

In our first paper (KSOl), we analyzed FUSE observations of the fluctuating He II 
absorption toward the bright (V = 16.1) quasar HE 2347-4342 at redshift z^^ = 2.885±0.005 
(Reimers et al. 1997; Smette et al. 2002). This object was discovered in the Hamburg/ESO 

survey for bright QSOs (Reimers, Kohler, & Wisotzki 1996) and is one of the brightest targets 
for FUV studies of He II, with specific flux Fa > 3 x 10"^^ erg cm^^ A^^ at 1180 A. The 
FUSE data include 376 ksec of observations between 1000-1187 A, taken at orbital night 
with a resolving power R = (A/AA) ^ 20,000. Using the H I Lya data from Kcck/HIRES, 
we analyzed both He II and H I Lya forests over the interval 2.3 < z < 2.9. About 50% of 
the He II features had H I counterparts with column densities Nhi > 10^^-^ cm~^. The He II 
to H I column density ratio rj ranged from 1-1000, with an average (77) f« 80. 

KSOl highlighted two major new results from the FUSE (He II) data: (1) the epoch of 
He II reionization appears to occur in this sightline at -ZHeii = 2.8 ± 0.2; (2) the range in the 
He II/H I ratio (rj) suggests that the ionizing background is highly variable throughout the 
IGM. Recent HST observations (Teller et al. 2002) show a broad distribution in QSO spectral 
indices, from — 3. Over the observed redshift interval, 2.3 < 2; < 2.7, the IGM toward 
HE 2347-4342 is "translucent" in He 11, with optical depths fluctuating between Tueii ~ 0.1 
and THeii > 3. Above z > 2.72, the IGM becomes relatively opaque in He II, except for 
three narrow transmission windows around 1160 A {z = 2.82) and 1175 A {z = 2.87) noted 
previously (Smette et al. 2002; KSOl). Below z < 2.72, the general trend is that the IGM 
becomes more transparent, with (thcIi) ~ 1, slowly decreasing at lower redshift. 

Except for the high-opacity regions {z = 2.7 — 2.9), the observed He II optical depths 
and He II/H I ratios are in reasonable agreement with models of IGM photoionization by a 
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QSO-type power-law spectrum, oc i/""", with (ag) = 1.8 ± 0.2 (Zheng et al. 1997; Telfer 
et al. 2002). Models of the He II/H I ratios in the IGM produced by QSO sources, filtered 
and reprocessed by the IGM, have been presented by FGS98 and Madau, Haardt, & Rees 
(1999). The high values [r] > 200) suggest additional contributions from softer ionizing 
sources (starburst galaxies) or heavily-filtered quasar radiation. 

In this second paper on HE 2347-4342, we provide an in-depth analysis and scientific 
discussion of these same FUSE He II data, in combination with new optical H I data taken 
with the VLT/UVES spectrograph (D'Odorico et al. 2000). We describe here the character- 
ization and subtraction of the background in the LiF (1000-1187 A) channels of FUSE. An 
independent analysis of both the LiF and SiC channels (920-1000 A) is described in a sep- 
arate paper (Zheng et al. 2003). We compare FUSE and HST/STIS data over the common 
interval from 1150-1187 A. We then analyze the He II and H I Lya absorbers, characterize 
their distributions in redshift and column density, and derive optical depths, THeii and tri. 
The observed absorption ratios, r) ~ 4rHeii/THi, can be used to understand the photoionizing 
spectra of the sources of the metagalactic radiation. We also explore the nature of the "void 
He II absorption", for which He II is an especially good probe. Finally, we address some of 
the open issues about the high-z IGM, including the He II reionization epoch, the ionization 
history of the IGM, and the range of ionizing spectra in the gaseous voids and filaments. 

The remainder of this paper is organized as follows. In § 2, we give details of our new 
FUSE data extraction and background subtraction. With the VLT data, we then derive the 
relevant parameters that characterize the two Lya forests (thi, thcii, i]). In § 3 wc discuss the 
scientific implications of these observations for the IGM and its ionization state a.t z < 2.9, 
just after the apparent He II reionization epoch. We conclude in § 4 with a summary of the 
highhghts of this new work and suggestions for future studies of He II with FUSE and HST. 

2. DATA ANALYSIS 

2.1. Provenance of the data 

As described in KSOl, HE 2347-4342 was observed with FUSE in two separate cam- 
paigns. The first observation, from 17 to 27 August 2000, comprised 351,672 s, 192,610 
s during orbital night. The second campaign, running from 11 to 21 October 2000, ac- 
cumulated 249,717 s, with 183,630 s during orbital night. For each observing campaign, 
HE 2347-4342 was centered in the 30" x 30" apertures. To maintain optical alignment of 
the four channels during these campaigns, FUSE was offset every other day to the nearby 
UV-bright star WD233 1-475 to perform routine adjustments to the mirror positions. FUSE 
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has an effective resolution R ^ 20,000 (Moos et al. 2000), which corresponds to a wave- 
length resolution element AA = 0.05 A at 1000 A. To avoid uncertainties in the effective 
spectrograph resolution and because the resolution element is highly oversampled, we bin 
our data to at least 0.05 A. This is approximately the minimum linewidth of the H I features 
we examine. 

We chose to neglect the wavelength range provided by the FUSE SiC detectors for 
two reasons. First, the SiC channels have a throughput that is ~3 times less than that of 
the LiF channels (A > 1000 A). Second, the higher effective area SiC channels (la and lb) 
appear near the bottom edge of the detectors, so that it is impossible to obtain a background 
spectrum based on unilluminated regions that bracket the spectrum. Recently, Zheng ct al. 
(2003) have succeeded in extracting a useful spectrum from the SiC channels by performing 
careful burst removal, pulse-height screening, and background modeling. 

HE 2347-4342 was observed during four nights (October 7-9 and October 19, 2000) using 
the UVES installed at the second VLT Unit Telescope (Kueyen). Twelve exposures were 
made in the dichroic mode using standard settings for central wavelengths of 3460/4370 A 
in the blue and 5800/8600 A in the red. The integration time for each wavelength range 
was 6 X 60 min. With a sht width of 1 arcsec, the resulting spectral resolution is about 
44,000. The raw frames were reduced at Quality Control Garching, using the UVES pipeline 
data reduction software (Ballester ct al. 2000). Finally, the individual vacuum-barycentric 
corrected spectra were coaddcd, resulting in an effective S/N ratio of roughly 100 per 0.05 A 
bin. When compared to the FUSE He 11 data, the H 1 spectrum is resamplcd to the same 
redshift vector and bin size; these sizes vary in our analysis from 0.05 to 0.2 A (see § 2.3). 

2.2. Data Calibration and Background Subtraction 

The faint flux levels of our target present unusual challenges to the instrumental cali- 
bration of FUSE. We can simply apply stock flux corrections and wavelength solutions to 
our data with high confidence that these solutions are stable in time and insensitive to the 
source fiux. However, the detector background correction that is usually applied to FUSE 
data by the CALFUSE pipeline (up to the latest version 2.2.1 as of August 2002) is in- 
sufficient to model correctly the spatial variations in the background for a source flux at 
(1 — 3) X 10~^^ erg cm~^ s" -1 A-\ In our 376 ksec of data, we collected roughly 5-20 source 
counts per 0.0067 A column in the extracted source spectrum, spread over 50 - 100 rows in 
the 2D detector image. This faint source level distributed over a large detector area leaves 
1-2 source -|- background counts per element in the illuminated regions of the detector image. 
This small source flux makes the background subtraction an unusually crucial piece of the 
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overall calibration. 

At the faint flux levels presented by HE 2347-4342, an accurate extraction of a background- 
subtracted, calibrated source spectrum from the two-dimensional recorded data required 
customized processing. We began our analysis with time-tagged photon lists representing 
all the data. We use only the nighttime data to avoid difficult problems related to scattered 
airglow light on the detectors. Up to the background subtraction, our analysis used the 
standard CALFUSE pipeline (v2.2.1) to process the photon lists, screen for and remove de- 
tector bursts, and create a 2D detector image. Wc applied the default CALFUSE day/night 
screening and truncated the range of detector pulse heights from channels 4 - 16 to minimize 
backgrounds. The resulting detector image has 16,384 "columns" (the x direction) and 1024 
"rows" (the y direction). These elements represent the digital quantization of the analog 
positions of photons striking the photocathode in the FUSE microchannel plate detectors 
(Sahnow et al. 2000). In this analysis, we avoid the term "pixels" because FUSE quantizes 
an analog quantity instead of using literal "pixels" as in CCDs. At this step, we depart from 
the CALFUSE procedure to produce a customized background subtraction with the scheme 
described below. 

Prior to constructing background vectors, for each channel/segment pair (LiFla, LiFlb, 
LiF2a, LiF2b), we added all the columns by summing in x. This summation reveals the y 
positions of the source-illuminated regions, including residual airglow lines present during 
orbital night. It establishes the height of the rectangular extraction windows and the y 
limits of the unilluminated background regions (see Table 1). Wc derive an x-dcpcndcnt 
background correction by taking unilluminated regions of the detector, above and below the 
illuminated source plus background region, and extracting a ID background spectrum as a 
function of x. This scheme gives us detailed information about the large spatial x variations 
in the background, but requires averaging over the variations in y that may occur beneath 
the illuminated portion. 

To construct the ID background spectrum for an extraction window of height n elements, 
we randomly drew n elements (at each x) from the background windows and summed over 
the number of elements in the extraction window. This sampling is done 300 times, and the 
final background spectrum is the average of the 300 trials at each x position. For example, 
for LiFla we randomly drew 80 bins from the 150 bins in the three available background 
windows. This sampling is done 300 times and the final background and uncertainty are 
drawn from the average and la variation in the distribution of samples. The final back- 
ground vector approximates the number of background counts collected in each row over 
the height of the extraction window. The random sampling accounts for fluctuations with 
y and assigns an error bar to the background subtraction, which is propagated through our 
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analysis. The final calibration steps apply the CALFUSE wavelength and flux solutions. 

For the FUSE data, we adopt the continuum determined from STIS/G140L data by 
KSOl, adjusted to account for the contributions (at the 5% level) of an intervening Lyman- 
hmit system at z = 2.739 with N(H I) = 2.5 x 10^'^ cm'^. This LL system makes a 5% 
contribution to the absorption at 1100 A. The extrapolation gives close agreement in He II 
transmission in the region (1157-1180 A) where our data and the STIS results (Smcttc ct 
al. 2002) overlap. The full recovery of He II to this continuum in transmission windows near 
1000 A gives us confidence that this extrapolation is acceptable at these wavelengths. 

The observations reported here push the limits of FUSE's design capabilities and exceed 
the limitations of the CALFUSE automatic data-reduction software. Our custom reduction 
was made as simple and accurate as possible, with careful attention to systematic errors and 
stability of the background subtraction. It appears that the differences between the KSOl 
reduction and the current work are comparable to the differences between overlapping seg- 
ments in our new extraction, which suggests that possible systematic errors associated with 
our background subtraction method and choice of extraction windows are well controlled. 
Furthermore, the differences between the October 2000 and August 2000 data and between 
LiFl and LiF2 arc consistent within statistical errors, and are not sensitive to the choices 
made in constructing the background spectrum. Table 2 shows the mean He II transmis- 
sion over several wavelength bands, the results of a cross-check with earlier studies by STIS 
(Smette et al. 2002) and FUSE (KSOl). 

Figure 1 shows an overlay of the FUSE (He II) data with the VLT/UVES (H I) data. 
The H I wavelength scale has been divided by a factor of 4.00178 (the ratio of He II and 
H I rydberg constants is not precisely 4.0) to align the H I and He II absorbers. One clearly 
sees the much larger He II absorption [r] ^ 1) predicted by theory. Other notable features 
include the transmission windows around 1132 A, 1160 A, and 1175 A, in close proximity to 
regions of high He II opacity. We also note the frequent voids in H I absorption. 

Figure 2 shows the broad-band transmission in the He II Lya forest from 1000-1180 A. 
As noted by KSOl, this shows the expected increase in broadband transmission towards 
shorter wavelengths. However, on a finegrained scale, AA = 0.05 - 0.20 A, the He II 
opacity shows substantial variations. In our analysis, we use various bin sizes to explore the 
dependence of H I and He II opacity on resolution. The smallest bin size we use is 0.05 A, 
or Az = 1.6 X 10"^, which corresponds to the resolution element of FUSE. 
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2.3. The He II/H I Ratio (ri) 



As defined, rj expresses the ratio of He II to H I column densities and can be used to infer 
properties of the ionizing radiation field. Owing to photon statistics and backgrounds, we 
believe that we can measure optical depths over the ranges 0.1 < THeii < 2.3 in the FUSE data 
and 0.02 < thi < 3.91 with the VLT/UVES. These values correspond to 10% uncertainties 
(FUSE) and 2% uncertainties (VLT/UVES) in measuring the optical depths of either strong 
or weak fines. In this paper, unlike KSOl, we do not make the assumption that He II 
absorption occurs in distinct lines. Instead of finding r) =N(He II)/N(H I) from the column 
density ratio, we use the observed opacities per resolution element to define r] — 4THeii/THi 
(the extra factor of 4 arises from bandwidth). As a consequence of these ranges for measuring 
optical depths, we can reliably determine values of r] between approximately 0.1 and 460. 

Figure 3 shows the observed distribution of r] with redshift, from z = 2.30 — 2.90, for 
two different wavelength bin sizes (0.05 A and 0.20 A). We also distinguish the absorbers 
in "voids" (thi < 0.05) from those in the H I "filaments" (stronger absorbers with Thi > 0.05). 
Note that rj is systematically higher in the voids. Binning the data from 0.05 A to 0.20 A 
reduces the number of points in the low (r/ = 0.1 — 10) and high (r] = 200 — 400) ends of the 
range, but it does not change the overall distribution. 

In practice, we must censor the data to avoid unphysical values of rj. First, we exclude 
points with Fx < 0. We also excise points within la of either the zero level or the continuum, 
on the grounds that they are indistinguishable statistically from and 1, respectively. Thus, 
the r] distributions shown below represent only a subset of the complete dataset. With these 
criteria for censorship, we find (77) 45. If the FUSE points within 1 cr of zero are used to 
provide lower limits on the He II optical depth, and thus lower limits on rj, a Kaplan-Meier 
survival statistics calculation implies (77) ~ 110. However, this neglects the upper limits 
available for 77. 

The uncensored mean value, (77) = 80, is in good agreement with theoretical predictions 
(FGS98) for an ionizing background produced by QSOs with mean spectral index (ag) — 1.8 
(Zheng et al. 1997; Teller et al. 2002). The wide range in 77, as well as the small-scale 
variations in redshift, were attributed in KSOl to the possible contributions of hard and soft 
sources to the ionizing background. However, this assertion deserves further scrutiny with a 
physical model of such sources and their spatial frequency. 

In photoionization equilibrium, the He II/H I ratio can be expressed (FGS98) as: 



77 = 



?^HeIII «HeII ^HI 
^HII ai^^ THell 



(1.70) 



Jhi (3 + ^4) 
Jueii (3 + ai) 



■T 



i0.055 
4.3 



(1) 
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Here, a^i , cn^Jn, Thi, and FHeii are the case-A recombination rate coefficients and pho- 
toionization rates for H I and He II, and Jhi and Jneii are the specific intensities of the 
radiation field at the H I (912 A) and He II (228 A) edges. The parameters ai and a4 
are the local spectral indices of the ionizing background at 1 and 4 ryd, respectively, which 
provide minor corrections to the photoionization rates. We adopt case-A hydrogenic recom- 
bination rates to H I and He II, appropriate for the Lya forest absorbers with very low 
neutral fractions. Over the temperature range 16, 000 K < T < 32, 000 K, we approximate 
«Hen = (1-36 X 10-12 ^^3 s-^)T^^-^'"' and ajjf = (2.51 x lO-^^ ^^3 s-^)T^^-^^\ where 
T4.3 = (T/10^-3 K). We also assume that nue/nn = 0.0789 {Y = 0.24) and that H and He 
are predominantly fully ionized (rie = l.lGriif). 

The metagalactic ionizing background, Ju{z), can be modeled (ShuU et al. 1999) from 
the intrinsic spectra of the ionizing sources (spectral index as) filtered through the IGM. 
Cosmological radiative transfer softens the spectrum and produces a steeper effective index, 
OefF > Q^s- For example, one can write the ratio of specific intensities at 1 and 4 ryd as 
Jm/Jneii = 4"°'', where a^s ^ 2.8 for = 2.0 (see Fig. 10 of FGS98). From equation (1), 
we then find (rj) ~ 80 for T^ ^ 1 and ai ^ a^. 

Figure 4 shows the observed range in ionizing spectral indices, ckeuV) oi radio-quiet 
QSOs from the HST spectral survey (Teller et al. 2002). This survey favors blue quasars 
with fairly clean lines of sight. This might bias the mean of the oeuv distribution to be 
bluer, and select against ones with very red cieuv- However, both radio-loud and radio- 
quiet samples show the same broad distribution of spectral indices. The upper curves in 
Fig. 4 show the spectral indices, as, required to reproduce the observed values of rj, using 
equation (1). We use bin- by-bin data for optical depths and convert them to the required a^, 
assuming 74.3 = 1 and ai — a^. In this approximation, Jm/^Heii = 4°", and 77 = 1.70 x 4"^^ 
from equation (1). It is also possible that some He II/H I variations could be produced 
by hot shocked gas (10^"^ K), for which 77 ranges from 10-500 in pure coUisional ionization 
equihbrium (Sutherland & Dopita 1993). 

We note that the range of ag (Fig. 4) is considerably larger than that of ckeuv, suggesting 
additional spectral softening. The absorbers in voids have higher spectral indices {{as) ~ 3.0) 
than those in filaments ((ctg) ~ 1.2), suggesting that the void absorbers may experience an 
ionizing radiation field softer than that of the "bare QSOs" . However, since r) has fine-grained 
variation over Az ~ 10~^, it may be inappropriate to use a mean metagalactic radiation field 
for all absorbers in the IGM. The ratio, Jhi/^HcII, must fluctuate due to "local effects", and 
the observed wide range of QSO spectral indices is proposed to account for most of these 
variations. This poses a paradox: the space density of AGN at z 2.5 is too sparse (mean 
distance of 30 Mpc between QSOs) to produce the observed 77-fluctuations {Az 10"^ or 
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l.ShjQ Mpc at z 2.5). We return to this issue in § 2.4. 

2.4. Statistics and Distribution of He II/H I Ratios 

In this subsection, we explore the variations of rj in greater detail, using an analysis of 
their correlations with voids and filaments in the Lya forest. The large number of 77 < 30 
points is unexpected, as it implies that there is a substantial volume in the universe where 
the extragalactic radiation background has a very hard spectrum. With a series of Monte 
Carlo calculations, we now explore the degree to which the high fraction of these points 
is due to bias in the points chosen to calculate 77. We begin with the H I spectrum and 
generate a new set of flux points based on the uncertainties given in that spectrum. We 
then convert this new set of flux points into an H 1 optical depth, thi- With the assumption 
that the underlying value of rj is actually a constant rjc, we then generate a corresponding 
set of He II optical depths, convert these optical depths into a set of flux points, and vary 
these points according to the uncertainty vector in our FUSE spectrum. We analyze this 
simulated spectrum in the same way as outlined for the observed spectrum, excising points 
with Fx < 0, and within 1 cr of zero or the continuum. Finally, we calculate a new set of 
values for 77. 

After performing a set of 1000 simulations, we arrive at the following conclusions: (1) 
For any assumed values of rjc, there is a trend to produce a simulated set of points with 

(77) < 77c; (2) Nevertheless, even for a value of 7/c chosen so that (7/) is the observed censored 
value (speciflcally, rjc = 52 and {rj) ~ 45 as in § 2.3), the fraction of low- 77 flux points in the 
data is larger than that in the simulations. In the data, 55% of all points have 7/ < 30, and 
30% of all points have 77 < 10. The mean fractions of such points in the simulations are 27% 
and 14%, respectively. No simulation had as many low-77 points as the data. 

Figure 3 indicates another check of the reality of the low 7/ points. We have made 
the same observational cuts in Thi for the points in the upper left and lower left plots. In 
the lower left panel, where the bin size is quadrupled and the S/N approximately doubled, 
the distribution in rj is much less likely to be affected by statistical fluctuations in the 
measured flux. The persistence of numerous low-77 points is reassuring to the interpretation 
that these points are real. As we will show, we believe the flne-scaled 77 variations are 
statistically signiflcant and arise from a combination of intrinsic variation in AGN spectral 
indices, coupled with radiative transfer effects in a non-uniform IGM. 

Although our analysis assumes that all absorption in the VLT spectrum is due to H I, 
and that all absorption in the FUSE spectrum is due to He II, both spectra are subject to 
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metal-line contamination. Songaila (1998) identified metal-line systems (C IV and Si IV) 
a.t z = 2.43814 and z = 2.63449, as well as metal lines connected with a Lyman-limit (LL) 
system (Reimers et al. 1997; Smette et al. 2002). The LL system has hydrogen continuum 
optical depth Tcont = 1-6 at Zi^i^ = 2.739, corresponding to N(H 1) = 2.5 x 10^'' cm"^ and can 
be seen in Figure 1 as the two saturated H 1 absorbers labeled on the He 11 scale at 1135 
and 1136 A (the corresponding He II lines fall within a detector gap, and so are not visible). 
In addition, the lack of a proximity effect in the He II Lya spectrum toward HE 2347-4342 
imphes a column density N(He II) = 2.4 x 10^^ cm~^ at Zabs ~ Zem (Reimers et al. 1997; 
Smette et al. 2002) [see § 3], which has highly ionized metal lines associated with it. Smette 
et al. (2002) report components with log N(C IV) = 12.8 and 12.1 as well as log N(0 VI) = 
13.6 associated with this absorption system. In general, for a range of assumed properties 
of this gas, the possible corresponding transitions are likely to be quite weak in the FUSE 
or VLT spectra. 

There also is the possibility that metal lines (C IV, Si HI) associated with lower redshift 
Lja absorbers may be present in portions of the VLT spectra. This is unlikely to impact our 
statistical conclusions strongly. At higher redshift, on average, (N(C 1V)/N(H I)) = 3 x 10^^ 
for N(H 1) ^ 10^^'^ cm~^. This would imply an equivalent width of only 1.2 mA for N(H I) 
pa 10^'^ '' cm~^, and much less for the more common lower column density absorbers. Finally, 
there is an unknown amount of H2 Lyman and Werner band absorption (A < 1120 A). 
This is a sightline with relatively low extinction, so that heavy Galactic H2 contamination 
is unlikely. We have explored the extent to which metal-line contamination would bias our 
statistical results. Simulating the effects of contamination, we randomize 10% and 20% of 
the opacity bins and reanalyze the data. In 100 realizations of each process, we have found 
that our statistical conclusions are robust. 

Next, we have examined the distribution of 77, both in regions of high H I opacity 
(denoted "filaments") and in low-opacity regions (denoted "voids"). Figure 5 shows the 
cumulative distributions of 1] in five different regions, two voids and three filaments, all in 
the redshift range 2.35 < z < 2.44. One can immediately see that the voids contain many 
more absorbers with large values of rj (soft ionizing sources) than the filaments. This effect 
was also noted in Figure 4. The radiation field in the voids could be dominated by starburst 
galaxies, whose massive stars lack the strong He II continua of QSOs. Alternatively, the 
AGN radiation reaching the voids may have been more strongly filtered by passage through 
dense filaments in the "cosmic web" . 

Figure 6 shows the distribution of optical depths, rneii, in two redshift intervals. In 
each panel, we show distributions for "voids" (thi < 0.05) and "filaments" (thi > 0.05) in 
the H I LycK forest. The fiat distributions in this logarithmic plot between 0.1 < Tneii < 1 
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are consistent with A'"(THeii) oc Tjjgjj, steepening at THeii > 1- This distribution is somewhat 
flatter than the corresponding H I column density distribution, /(A^m) oc N^i, with /3 = 
1.4 - 1.5 for 1.5 < z < 4 (Kim et al. 2001). 

Figure 7 shows a closeup of one of the "void" regions {z — 2.40 — 2.42) in which the He II 
transmission is unusually large, again with many high-77 absorbers suggestive of local soft 
sources. However, the difficulty with this simple explanation is that the expected high space 
density of (dwarf starburst) galaxy sources should produce a uniform ionizing background, 
not the observed fluctuating fleld of rj values. An alternative interpretation is that the hard- 
spectrum QSOs might preferentially be located in the overdense regions associated with 
H I filaments. In this case, the filament absorbers would see a harder spectrum (lower a^, 
lower 77) than the voids. Ultimately, these speculations need to be confirmed by cosmological 
radiative transfer calculations in an inhomogeneous IGM. 



3. DISCUSSION 

In this paper, we have used FUSE (He II) and VLT/UVES (H I) data to analyze 
the fiuctuating intergalactic absorption between z = 2.3 — 2.9 toward the high-redshift QSO 
HE 2347-4342. This interval appears to span the He II reionization epoch at ^Heu = 2.8±0.2, 
which is considerably lower than that, zm — 6.2 ± 0.2, estimated for hydrogen reionization. 
For a standard WMAP cosmology (Spergel et al. 2003), with ifo = 71 ± 4 km Mpc-^ 
Qrn — 0.27 ± 0.02, and I^a = 0.73 ± 0.02, the interval of time between these redshifts 
corresponds to At — (1.45 ± 0.13)h^Q Gyr. 

As in KSOl, we confirm the slow reduction of THeu at lower redshifts, z < 2.7, in 
good agreement with theoretical expectations (FGS98) of an IGM photoionized by QSOs. 
The He II "reionization epoch" in this one sightline is characterized by several regions of 
high optical depth between z ^ 2.75 — 2.85, separated by narrow transmission windows at 
z = 2.82 and 2.87. Clearly, FUSE and/or HST studies of additional sightlines are needed to 
understand the cosmic variance in these results. 

Our major findings on IGM structure and sources of ionizing intergalactic radiation are: 

1. We see fine-grained variations in He II/H I ratios (77) on scales Az ~ 10~^, correspond- 
ing to physical lengths ~ 1.3h^Q Mpc at ^ a; 2.5. 

2. The 77- variations appear to refiect "local radiation effects" from QSOs whose spectral 
indices range from as — — 3. A few regions with 77 < 30 could lie close to hard 
ionizing sources. 
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3. The He II/H I ratio appears to be larger in voids in the H I Lya forest {thi < 0.05) 
compared to filaments. In these voids, up to 80% of the pathlength has rj > 100. 

These results suggest that, if both He 11 and H 1 are photoionizcd, the mctagalactic radiation 
field is highly variable and strongly affected by local sources. In principle, the ?7-fluctuations 
could arise solely from variations in IGM opacity. But, in view of the broad observed range 
in AGN spectral indices, ctg = — 3, it seems probable that the fiuctuations arise from 
spatial variations in the ionizing fiuxes ("local effects"). Moreover, observations of the Lya 
absorbers at z = 1 — 3 (Kim et al. 2001) suggest a low spatial frequency of strong absorbers 
that dominate the opacity at the H 1 or He 11 ionization edges. Optical depths unity in these 
continua correspond to H 1 column densities lO^^'^ cm~^ and 10^^'^(100/?7) cm~^, respectively. 
Such systems are too rare in the Lya forest at z ^ 2.5 to produce the fine-grained variations 
over Az pa 10"^. 

However, the simple explanation that QSO emissivity variations translate into fluctua- 
tions in rj needs a more flrm physical basis. The variable-index QSOs at z = 2.3 — 2.8 must 
produce spatial variations in the radiation field over Mpc-scales {Az = 10^'^), whereas the 
co-moving space density of bright (L*) QSOs at z = 2.5 is 

0QSO ~ 10"^ h^o Mpc"^ which 
corresponds to mean physical separation, ((igso) = 0qso/(1 + ^) ~ '^Ohj^ Mpc. As discussed 
by Giroux, Fardal, & ShuU (1995), the "proximity spheres" for bright QSOs are ~1 Mpc 
at 1 ryd but as high as 35 Mpc at 4 ryd. At these distances, the QSOs's local radiation 
field matches the metagalactic background and can create detectable changes in the H I and 
He II line opacity. In the case of HE 2347-4342, Smette et al. (2002) noted the absence of 
any observed proximity effect in He 11 absorbers near Zcm = 2.885. They estimated that 
a column density N(He 11) > 2.4 x 10^^ cm^^ was required to absorb the QSO's ionizing 
radiation at 4 ryd. Thus, the QSO radiation field in the He 11 continuum may be strongly 
affected by absorption from both the galactic nucleus and the IGM. 

If a significant fraction of AGN exhibit He 11 intrinsic absorption, then the ratio of 
fiuxes at 1 and 4 ryd can be larger than implied by the range in q;euv (Telfer et al. 2002), 
which is only fitted between 1200 and 500 A. Moreover, if this intrinsic absorption partially 
obscures the nuclear continuum in only some directions, other QSOs off the sighthne could 
create regions of hard and soft ionizing radiation. These AGN could produce fine-grained 
variations in rj, owing to their own patchy absorption. This idea was proposed for a "helium- 
reionizing quasar" 3.2 Mpc off the sighthne to Q0302-003 (Jakobsen et al. 2003). 

We now can characterize the difficulty with the simple explanation of attributing fluc- 
tuations in rj to variable QSO spectral indices. The FUSE/VLT data suggest that the local 
radiation held varies on 1 Mpc scales, whereas bright QSOs are 30 Mpc apart, on average. 
However, their proximity spheres extend to 1 Mpc at 1 ryd, and to 35 Mpc at 4 ryd. The 



-14- 



observed fluctuations in rj are produced by the combination of the wide range in source 
emissivities and AGN/IGM opacities at 1 and 4 ryd. The non- uniformity of the metagalac- 
tic background, J,^, needs to be confirmed by numerical simulations that carefully treat the 
radiative transfer (absorption, filtering, reprocessing) of ionizing radiation through the cos- 
mic web and include the range of QSO spectral indices. In standard models of the ionizing 
background (FGS98; Madau et al. 1999), the radiation field at a given point is produced 
by many ionizing sources within a sphere of radius Vq, the comoving attenuation length for 
"optical depth unity" in the continuum. Fluctuations in the ionizing background are then 
set by statistics of sources within that sphere. 

Fardal k ShuU (1993) estimated that iYgso = 300 - 1000 QSOs lie within this "attenu- 
ation sphere" at z fa 2.5 for a standard QSO luminosity function (Boyle et al. 1991). From 
these large numbers, one naively expects the H I ionizing background to be fairly uniform. 
However, these estimates differ at 1 ryd and 4 ryd. FGS98 found that, at z 2.5, the H I 
and He II continuum opacity gradients, dTcont/dz are 2.5 (H I) and 10 (He II). Therefore, 
the spheres of optical depth unity extend over Az fa 0.4 (H I physical length 140 ± 10 Mpc) 
and a; 0.1 (He II physical length 34 ± 1 Mpc). 

These results predict much larger intensity fluctuations in the He II continuum than in 

the H I continuum, both because the He II attenuation sphere is smaller, and because the flat- 
spectrum QSOs (fts = — 1.5) are rarer than the the steep-spectrum QSOs (a^ = 1.5 — 3.0), 
an effect seen clearly in Figure 4 and in Tclfcr et al. (2002). In numerical hydrodynamic 
simulations, the IGM is distributed non-uniformly in filaments, and the spherical approxi- 
mations break down. The propagation of 1 ryd and 4 ryd radiation differ considerably, owing 
to the frequency dependence of the opacity. More sophisticated models should be able 
to confirm our hypothesis that the observed 77-fluctuations are the result of non-uniformities 
in opacity, coupled with the observed range in QSO spectral indices. 



4. SUMMARY OF SCIENTIFIC IMPLICATIONS 

This paper continued the study of the He II/H I absorption toward HE 2347-4342, 
using both FUSE and VLT spectra to derive values of ?] = 4rHeii/'rHi ranging from 0.1 to 
460. The combination of He II and H I absorption data allows us to probe the metagalactic 
ionizing spectrum and its sources. The large observed fluctuations in the H I and He II Lya 
absorption, and the wide range (0.1-460) in the ratio 77 = N(He II)/N(H I) suggest a wide 
range of ionizing sources and filtering of the ionizing spectra by a filamentary IGM. The fine- 
grained variations in rj indicate that the ionizing background is far from uniform, probably 
due to variations in both source emissivity (spectral index) and opacity fiuctuations. A 
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recent HST survey (Telfer et al. 2002) of QSO ionizing spectra suggests that most of the 
observed variations (20 < r] < 200) can be understood through the wide range in QSO 
spectral indices (0 < ct^ < 3). 

The most significant scientific results of our study are: 

• The observed fine-scaled {Az = 10~^) variations in r], on physical scales ~ l.Sh^Q Mpc 
are much smaller than typical distances between AGN or Lyman-limit absorbers. 

• The wide range in rj (from 0.1-460) implies a broad distribution of effective ionizing 
spectral index, a^, ranging between —1 and 4 (F,^ oc i/~"^). This range exceeds the 
observed range of AGN spectral indices, suggesting additional spectral softening due 
to radiative transfer through the AGN and a filamentary cosmic web. 

• An observed systematic correlation of higher r] (softer radiation fields) for absorbers in 
voids in the H I distribution (thi < 0.05). 

• An inferred distribution in He II line optical depth, N(THeii) oc r^^n between < 

THell < 1- 

We conclude with a discussion of future needs and challenges. The data on HE 2347-4342 
(-^em = 2.885) provide important evidence on the epoch of He II reionization {zr — 2.8 ±0.2), 
and on the spectral variations and spatial distribution of ionizing sources. Nevertheless, 
one always worries whether cosmic variance might produce differences along other sight- 
lines. It would be valuable to obtain He II absorption data along several other FUSE or 
HST sightlines, over the critical redshift range 2.4 < z < 3.2. At least one other sightline 
(HS 1700+6416 at Zem = 2.72) has just been observed with FUSE, and the results look 
promising. 

One of the major theoretical challenges is to understand why the He II reionization 
epoch (zhcII = 2.8 ±0.2) occurs much later than the H I reionization epoch (^hi = 6.2 ±0.2). 
Did He II experience "double re- ionization" , after an initial ionization epoch at ^ = 6 — 10 
by Pop HI, low-metallicity stars (Venkatesan, Tumlinson, & ShuU 2003)? Vestiges of this 
effect might still be visible in low-density voids, although much of the effect will be washed 
out by the rapid increase in QSO fluxes at 2; < 6. The observed high rj that we find in the 
voids would seem to rule out local hard spectra and relic ionization from early metal-free 
star formation. 

Our finding that 77 varies on small (1 Mpc) distance scales poses a quandary for the 
nature of the ionizing sources. Although it is tempting to conclude (KSOl) that the wide 
range of AGN spectral indices can explain much of the He II/H I variation (20 < 77 < 200), 
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one still must confirm, through numerical simulations of ionizing sources embedded in the 
cosmic web, that these 77-variations actually occur naturally. 

Another intriguing result from our study is that 77 is systematically larger in voids 
(thi < 0.05), compared to filaments of H I Lya absorbers. Does a higher 77 really imply a 
softer ionizing radiation field, characteristic perhaps of starburst galaxies? Can we therefore 
use He II absorbers to probe the distribution of gas in voids? 

If the variations rj can be explained by local fluctuations in the ionizing radiation field, 
enhanced by radiation transfer through the IGM, there may be additional signatures that 
take place around the epoch of He II recombination, ^Heii = 2.8 ± 0.2. It is still not certain 
that the effect is entirely the result of variable-index QSOs; soft spectra can also arise from 
starburst galaxies. Deep imaging of these fields at 2.4 < z < 2.8 might detect the actual 
sources of ionizing radiation. At z ~ 2.5, an L* galaxy (m^ = —19.5, with 21 Gpc luminosity 
distance) would have magnitude ~ 27. One would need to search angular scales out to 
2' from the sighthne (1 Mpc offset at z = 2.5). One might also hope to obtain additional 
physical information on the ionization state of these H I and He II absorbers by looking for 
the corresponding C IV, Si IV, and O VI absorption in the NUV and optical. 
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Johns Hopkins University and from NASA grant NAG5-13013 at the University of Colorado. 
J. M.S. and M.L.G. acknowledge support at the Colorado astrophysical theory program from 
NASA LTSA grant NAG5-7262 and NSF grant AST02-06042. 



Table 1. Background Extraction Windows'^ 



Segment Extraction Window Background Windows 



LiFla 


530 - 


609 


450 - 499, 720 - 769, 860 - 909 


LiFlb 


470- 


607 


300 - 449, 650 - 799 


LiF2a 


660 - 


709 


600 - 659 


LiF2b 


690 - 


729 


600 - 636, 734 - 770 



^We sum data from all columns (x) between the listed columns 
(y values above). 
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Table 2. Measured Transmission in He II 



Band (A) 


Kriss et al. (2001) 


Smette et al. (2002) 


This work 


1157.75-1158.25 


0.69 ± 0.06 


0.78 ± 0.09 


0.80 ± 0.07 


1159.0-1161.0 


0.68 ± 0.03 


0.64 ± 0.04 


0.63 ± 0.04 


1161.5-1168.3 


0.02 ± 0.02 


0.06 ± 0.02 


0.04 ± 0.02 


1174.0-1175.0 


1.03 ± 0.06 


0.93 ± 0.04 


0.94 ± 0.06 


1175.1-1179.3 


0.07 ± 0.03 


0.02 ± 0.02 


0.00 ± 0.03 
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Fig. 1. — Overlay of the Lya absorption of He II (from FUSE - top) and of H I (from 
VLT/UVES - bottom). Wavelengths for the lower panels (H I data) are divided by 4.00178 
to align them with the He II data. The extrapolated STIS continuum is overlaid for the 
FUSE data. In the VLT data, "filament" points with thi > 0.05 are shown in red and 
blocked in yellow above. 
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Fig. 2. — Flux transmission through FUSE He II absorption in several defined windows. Note 
the three windows of high transmission between z = 2.7 — 2.8, fianked by low-transmission. 
The fiux has a slow recovery of transmission at ^ < 2.4. Table 2 gives a comparison with 
previous studies. 
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Fig. 3. — Left panels show the distribution with redshift of rj = N(He II)/N(H I), determined 
here &sr] = 4rHeii/THi, in wavelength bins of 0.05 A (top panels) and 0.20 A (bottom panels). 
The distribution of rj, integrated over 2.3 < z < 2.9 is shown in the right panels. With the 
accuracy of these data, we can reliably measure optical depths in the ranges 0.1 < thcii < 2.3 
and 0.02 < thi < 3.91. Gap at z = 2.6 lies between the LiFa and LiFlb FUSE detector 
segments. "Filaments" in the Lya forest are plotted in black (thi > 0.05), while "voids" 
have 0.02 < thi < 0.05 and are plotted in red, some as lower limits. The large fluctuations 
in T] suggest wide variation in the spectra of the ionizing sources, density fluctuations, or 
significant effects of radiative transfer in the IGM. 




2 



Fig. 4. — (Bottom curve) Observed distribution of the observed QSO ionizing spectral in- 
dices, oeuv, for radio-quiet AGN (Telfer et al. 2002). The other curves show the spectral 
indices, a^, required (using eq. [1]) to reproduce the observed values of He II/H I ratios. Note 
the systematically softer radiation fields (high 77, high a^) in the voids (Red: thi < 0.05) 
compared to filaments (Green: thi > 0.05). The breadth of these curves and the clear offset 
between voids and filaments suggest that additional spectral softening may be present. 
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Fig. 5. — Cumulative distribution functions of log 77. For regions with significant H I opacity, 
roughly 40% of the bins have log?] > 2.0. For the void region, 80% of the bins have 2.0 < 
log?7 < 3.0. Uniformly high 77 suggests the influence of local soft sources (e.g., widely 
distributed dwarf galaxies) or heavy filtering by the IGM. 
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Fig. 6. — Distribution of optical depths, rneii, in two redshift intervals, 2.29 < z < 2.56 
(regions of lower optical depth) and 2.62 < z < 2.88 (regions of patchy, variable opacity. In 
each panel, we show distributions for "voids" and "filaments" in the H I Lya forest. The 
distributions look flat (per logrneii, between 0.1 < rneii < 1, and steepen at THeii > 1- This 
corresponds to N{t) oc r^^ 
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Fig. 7. — Closeup of filaments and voids (z = 2.36 - 2.44). In the voids, (thi < 0.05), the 
He II transmission is large and rj is uniformly high, perhaps indicating the influence of local 
soft sources. 



